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Preface 


This  thesis  is  concerned  with  the  development  of  a computer 

code  to  analyze  alpha  particle  spectra.  The  reference  peak  used 
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in  the  program,  obtained  from  a Po  source,  was  expressed  as 
a table  of  values  specifying  the  peak  shape  channel  by  channel. 
The  results  of  the  analysis  were  not  within  the  limits  desired 
but  this  is  thought  to  be  a result  of  different  source  thick- 
nesses used. 
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Abstract 


An  existing  computer  code,  written  to  analyze  pulse-height 
spectra  from  charged  particles  measured  with  a mass  spectro- 
graph, was  modified  to  analyze  pulse-height  spectra  detected 
by  semiconductor  detectors.  The  program  fits  the  spectra  with 
a reference  peak  which  is  expressed  as  a table  of  values  speci- 
fying the  peak  shape,  channel  by  channel.  The  modification 
involved  increasing  the  size  of  the  reference  peak  to  20^.8 
channels  to  allow  analysis  of  pulse-height  spectra  from  alpha 
particles  with  energies  between  1;  and  6 MeV.  The  results  of 
the  analysis  were  strongly  dependent  upon  the  reference  peak 
chosen;  hence,  reference  peak  selection  became  the  main  emphasis 
of  this  report.  The  three  reference  peaks  used  were  obtained 
from  (1)  an  empirical  function,  (2)  a single  alpha  peak  obtained 
from  Po  , and  (3)  an  entire  spectrum  of  an  isotope  to  be 
analyzed.  It  was  desired  that  the  error  resulting  in  determining 
the  area  be  within  three  standard  deviations  of  the  actual  area. 
The  error  resulting  in  each  analysis  was  greater  than  the  pre- 
scribed value  and  was  suspected  to  have  occurred  because  of  mis- 
calculation of  the  reference  peak. 


vii 


A COMPUTER  CODE  TO  ANALYZE  ALPHA  SPECTRA 
I . Introduction 

This  report  presents  the  results  of  an  investigation  of 
methods  for  the  analysis  of  pulse-height  spectra  produced  by 
alpha  particles  detected  by  semiconductor  detectors.  Several 
goals  explored  activities  of  radionuclides,  resolution  of  spectra, 
the  effects  of  parameters  such  as  source  thickness  and  width 
of  pulse-height  interval  in  obtaining  the  pulse-height  spectra, 
and  various  types  of  algorithms  to  analyze  pulse-height  spectra. 
Major  emphasis  was  placed  on  obtaining  a peak  shape  represen- 
tative of  a single  alpha  peak  and  using  this  as  a reference  peak 
along  with  existing  algorithms  to  fit  alpha  spectra. 

Background 

The  general  technique  of  pulse-height  analysis  consists  of 
determining  the  area  under  each  peak  in  the  pulse-height  spec- 
trum. The  area  is  frequently  determined  by  obtaining  an  analyt- 
ical expression  which  will  reproduce  the  peak  shape.  The  ex- 
pression is  fit  to  the  data  by  the  method  of  least  squares;  then 
numerically  integrated  to  obtain  the  area.  The  area  represents 
the  total  number  of  particles,  with  a certain  energy,  detected 
by  the  equipment. 

The  method  of  pulse-height  analysis  has  been  applied  to 
gamma-ray  spectra  for  quite  some  time.  In  gamma-ray  spectra, 
the  photopeaks  are  approximated  by  Gaussian  functions,  which 
can  be  numerically  integrated  without  difficulty  to  provide 
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precise  determination  of  the  areas  under  the  photopeaks. 

Alpha  particles  unlike  gamma-rays  lose  energy  in  small 
decrements  as  they  collide  with  the  atoms  in  the  detector.  As 
a result  cf  the  possibility  of  small  energy  losses,  the  peaks 
in  the  pulse-height  spectrum  are  no  longer  characterized  by  a 
pure  Gaussian  function,  but  rather  by  a Gaussian  shaped  leading 
edge  and  an  exponential  shaped  trailing  edge.  It  is  difficult 
to  obtain  an  analytical  expression  for  the  resulting  peak  shape. 

Currently,  a crude  alpha  pulse- height  analysis  technique 
is  being  used  to  analyze  alpha  spectra  obtained  from  sources, 
such  as  Pu239  and  Pu2^-0.  The  transition  energies  of  the  alpha 
particles  emitted  by  these  isotopes  differ  by  only  a few  keV; 
therefore,  semiconductor  detectors  are  not  able  to  resolve  the 
peaks  resulting  from  each  isotope.  However,  the  contribution 
from  each  isotope  can  be  determined  by  using  the  sum  of  the 
counts  from  both  isotopes  along  with  the  known  relative  amounts 
of  each  isotope  obtained  by  mass  spectrometry.  The  limiting 
factor  in  this  procedure  is  the  determination  of  the  absolute 
sum  of  the  counts  from  both  isotopes.  A summation  region  ap- 
proach is  currently  being  used  to  determine  the  sum  of  the  counts. 
This  method  consists  of  determining  the  area  under  the  peak  by 
summing  the  counts  per  channel  for  all  channels  within  a spec- 
ified region  containing  the  peak.  A serious  drawback  in  this 
method  results  from  information  lost  in  the  tails  of  the  peaks. 
This  occurs  because  of  the  necessity  to  specify  a region  into 
which  the  peak  must  fall.  The  regions  chosen  are  mutually  ex- 
clusive; hence,  counts  found  in  the  long  tail  extending  beyond 
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the  specified  region  are  lost.  Greater  precision  could  be  ob- 
tained if  the  alpha  spectra  could  be  analyzed  using  a peak  shape 
which  would  closely  approximate  the  actual  alpha  peak  shape. 

Problem  Statement 

The  object  of  this  thesis  was  to  modify  an  existing  com- 
puter code  and  to  obtain  a peak  shape  representative  of  a single 
alpha  peak.  This  peak  shape  was  then  used  as  a reference  peak 
by  the  program  to  analyze  the  alpha  spectra. 

Scope 

The  study  was  restricted  to  the  use  of  an  existing  computer 
code  which  was  modified  to  allow  analysis  of  alpha  spectra.  A 
reference  peak  representative  of  the  peaks  to  be  analyzed  was 
required  by  the  program  to  perform  the  analysis;  hence,  a large 
amount  of  effort  was  devoted  to  obtain  a representative  alpha 
peak. 

General  Approach 

The  program  used  in  this  study  was  obtained  from  a program 
called  AUTOPIT,  written  by  J.  R.  Comfort  of  Argonne  National 
Laboratory.  The  major  modification  was  to  increase  the  allowed 
number  of  channels  in  the  reference  peak  from  IpO  channels  to 
2048  channels.  After  modification,  the  program  was  renamed 
ALPHAFIT  and  was  tested  to  determine  its  peak  fitting  capability. 

An  analytical  function  was  used  to  generate  emulated  pulse- 
height  spectra.  For  simplicity,  a Gaussian  pulse-height  dis- 
tribution was  used.  Studies  were  made  of  the  capabilities  of 


the  computer  program  to  resolve  two  peaks  whose  relative  ampli- 
tudes and  positions  were  varied.  Other  variables  studied  were 
the  effects  of  resolution  and  magnitude  of  pulse-height  inter- 
val used  to  generate  the  simulated  spectra.  Two  series  of  tests 
were  performed.  Statistical  fluctuations  of  the  individual 
points  in  the  simulated  spectra  were  included  in  the  second 
series  of  tests,  but  not  in  the  first  series  of  tests. 

Several  methods  were  used  to  obtain  a reference  peak  which 
would  be  representative  of  the  alpha  peaks  to  be  analyzed.  The 
first  method,  developed  by  P.  De  Regge  (Ref  2:271-272),  allowed 
a single  alpha  peak  to  be  extracted  from  a complex  alpha  spec- 
trum such  as  that  obtained  from  A m^1 . 

• The  second  method  involved  using  a single  alpha  peak  ob- 
21 0 

tained  from  Po  directly  as  a reference  peak. 

The  last  method  was  based  on  using  an  entire  spectrum  as 
a reference  peak.  Since  the  types  of  isotopes  in  the  spectra 
are  known,  a reference  peak  composed  of  an  entire  spectrum  of 
each  isotope  in  the  spectra  can  be  used  as  a reference  peak 
for  that  isotope.  The  resulting  analysis  should  be  more  precise 
and  less  time-consuming  than  using  single  alpha  peaks. 
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II.  Theory 


Before  becoming  involved  in  the  techniques  used  to  deter- 
mine the  area  under  the  peaks  in  a pulse-height  spectrum,  it 
is  important  to  know  the  cause  and  meaning  of  the  pulse-height 
distributions.  The  information  presented  in  this  section  de- 
scribes the  source  of  alpha  particles  and  the  way  they  interact 
with  the  detector  to  produce  the  pulse-height  distributions. 

A method  on  how  to  relate  the  area  under  the  pulse-height  spec- 
trum to  the  activity  of  the  source  is  described  followed  by  a 
description  of  the  algorithm  used  by  the  program  to  minimize 
the  error  between  the  reference  peak  and  the  actual  spectrum. 

Alpha  Decay  Process 

The  motivation  behind  the  development  of  a computer  code 

to  analyze  alpha  spectra  was  the  need  to  analyze  alpha  spectra 

produced  by  mixtures  of  low-activity  sources,  such  as  plutonium, 

americium,  and  curium.  However,  during  the  development  of  the 

program,  it  was  necessary  to  use  sources  of  sufficient  activity 

to  reduce  counting  times  and  to  increase  counting  statistics. 

The  two  main  sources  used  in  the  developmental  stage  were  Am^^ 
21 0 

and  Po  . These  isotopes  decay  by  alpha  particle  emission 
either  to  the  ground  states  or  the  excited  states  of  their 
daughter  nuclei.  Figure  1 shows  the  simplified  alpha  decay 
scheme,  describing  the  major  branchings  (Ref  4:435»U03)» 


Alpha  Particle  Interaction 

When  alpha  particles  pass  through  matter,  they  lose  en- 
ergy by  excitation  and  ionization  of  the  atoms  of  that  matter. 
The  mechanism  responsible  for  this  energy  loss  is  the  inter- 
action of  the  Coulomb  fields  of  the  particle  with  those  of  the 
bound  electrons  in  the  material.  Deflection  of  the  incoming 
alpha  particles  is  negligible  because  of  the  large  mass  of  the 
• alpha  particle  relative  to  the  mass  of  the  electrons. 

As  an  alpha  particle  passes  through  a medium,  the  energy 
dE  transferred  to  that  medium  by  ionization  and  excitation 
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along  a differential  path  length  dx  is  given  by  the  following 
theoretical  expression  (Ref  5*3-14-): 

-dE/dx  = kne^z2m/mv2  (l) 

where  B is  defined  as  the  atomic  stopping  number  and  is  given 
by 

B = Z[ln(2mv2/I)  - ln(1-p2)  - p2]  (2) 

with 

e = electronic  charge 
m * electron  rest  mass 
E = kinetic  energy  of  primary  particle 
ze  = charge  of  primary  particle 
v - velocity  of  primary  particle 

•i 

N * number  of  absorber  atoms  per  cmJ 

Z = atomic  number  of  absorber 

p * v/c,  where  c = velocity  of  light  in 
a vacuum 

I « geometric-mean  ionization  and  excitation 
potential  of  absorbing  atoms;  cannot  be 
calculated  accurately  and  generally 
regarded  as  a constant  for  each  element. 
The  value  of  I must  be  experimentally 
determined  for  each  element. 

In  the  nonrelativistic  energy  range,  the  rate  of  energy 
loss  -dE/dx  depends  on  the  velocity  of  the  alpha  particle 

p 

through  the  1/v  term  in  Eq(1).  The  increase  in  the  rate  of 
energy  loss  with  decrease  in  velocity  is  to  be  expected.  At 
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lower  velocity  the  alpha  particle  spends  more  time  in  the 
vicinity  of  bound  absorber  electrons,  thereby  increasing  the 
probability  of  ionization  and  excitation. 

In  the  relativistic  energy  range,  the  rate  of  energy 
loss  -dE/dx  passes  through  a minimum  which  is  followed  by  a 
slow  increase  with  increasing  particle  energy. 

The  ranges  of  alpha  particles,  which  are  initially  mono- 
energetic,  are  not  precisely  the  same.  This  variation,  or 
straggling,  is  due  to  the  statistical  nature  of  the  energy- 
loss  process.  The  energy  losses  which  bring  the  alpha  particles 
to  rest  consist  of  a large  number  of  individual  energy-transfer 


events  of  varying  magnitudes.  Most  of  the  energy  transfers 
are  relatively  small.  The  statistical  nature  of  the  energy- 
loss  process  allows  variations  in  the  number  of  events  oc- 
curring per  unit  path  length  and  variations  in  the  amount  of 
energy  transferred  per  event.  Pew  events  involve  large  amounts 
of  energy  transfer. 

The  absorption  of  alpha  particles  is  usually  studied 
experimentally  by  measuring  the  number  of  ion  pairs  produced 
per  unit  path  length  of  travel.  This  quantity  is  commonly 
referred  to  as  the  specific  ionization.  The  rate  of  energy 
loss  is  related  to  the  specific  ionization  through  a quantity 
w,  which  is  the  ratio  of  the  energy  lost  by  a charged  particle 
/ - to  the  total  ionization  produced  by  that  charged  particle. 

Values  of  w depend  upon  the  material  and  the  state  of  the 
material  with  which  interaction  occurs,  and  the  type  and  the 
energy  of  the  incident  particle. 
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Alpha  Pulse  Height  Distributions 

Alpha  particles  incident  on  solid  state  detectors  deposit 
energy  within  the  charge-depletion  region  of  the  detectors. 

The  energy  deposited  by  the  alpha  particles  produces  a number 
of  charge  carriers  (electron-hole  pairs)  proportional  to  the 
energy  deposited.  These  free  charge  carriers  which  are  swept 
from  the  depletion  region  produce  a pulse  height  which  is 
proportional  to  the  number  of  electron-hole  pairs  collected, 
and  therefore  proportional  to  the  energy  deposited  by  the 
incident  radiation.  The  pulse  height  out  of  the  detector  can 
be  expressed  by  the  following  equation: 


Ph  = [(Edep^w)  x e x y]  x 1/0  (volts)  (3) 

with 

E,  = energy  deposited  by  the  alpha 

p particle  (eV) 

w = average  energy  required  for  the 
production  of  an  electron-hole 
pair  (eV/e-h  pair) 

e = electron  charge  (coulombs) 

y = charge  collection  efficiency. 

This  is  essentially  1 in  detectors 
made  of  high  quality  single  crystals 

N = Edep^w  * num^er  electron-hole 
pairs  produced 

q = Ney  = total  charge  collected  (coulombs) 

C = capacitance  of  the  detector  (picofarads) 


-•  ■ :■ 
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Eq  (3)  can  be  rewritten  as 

ph  = q/C  (volts)  (4) 

the  output  of  the  linear  amplifier  is 

PH  * q x G (volts)  (5) 

q - total  charge  collected 
(coulombs ) 

G = gain  of  linear  amplifier 
(volts/coulomb) 

A multichannel  analyzer  is  used  to  count  the  number  of 
times  pulses  with  amplitudes  between  h and  h + Ah  occur  during 
the  acquisition  time.  This  information  is  displayed  on  ?.n 
oscilloscope  with  the  pulse-height  intervals  plotted  on  the 
abscissa  and  the  number  of  times  each  pulse-height  is  observed 
on  the  ordinate.  Since  the  energy-loss  process  is  statistical 

* 

in  nature,  the  resulting  pulse-heights  will  form  a pulse-height 
distribution.  The  energy  of  the  alpha  particles  are  associated 
with  the  location  of  the  maximum  in  the  pulse-height  distri- 
bution. The  distribution  of  pulse-heights  does  not  reflect  a 
distribution  of  energy  lost  by  the  alpha  particles  because  each 
particle  loses  an  energy  EQ  in  the  detector.  The  tailing  effect 
on  the  low  energy  side  of  the  pulse-height  distribution  occurs 
when  the  energy  of  the  alpha  particles  is  degraded  by  passage 
through  material  in  the  source  and  the  dead  layer  of  the  detector. 


The  pulse-height  at 
given  by 

with 


The  length  of  the  path  and  thus  the  loss  in  energy  by  each  alpha 
particle  depends  on  the  angle  at  which  it  is  emitted  from  the 
source  and  the  extent  to  which  it  deviates  from  normal  incidence 
when  it  passes  through  the  dead  layer  of  the  detector.  The 
effects  of  various  path  lengths  of  alpha  particles  through  gold 
can  be  seen  in  Figure  2. 


Figure  2.  The  oc-particle  Energy  Distribution  for  a 
Thin  Source  of  Pu  After  Transmission  Through 
Various  Thicknesses  of  Gold  (Ref  6:108) 


Electronic  noise  in  the  detection  system  will  further  in- 
crease the  spread  of  the  peak  in  the  pulse-height  spectrum. 
Electronic  noise  can  be  determined  by  measuring  the  spread  in 
the  pulse-height  spectrum  from  a pulser  signal  which  is  equal 
to  the  average  pulse-height  produced  at  the  input  to  the 


preamplifier  by  the  source.  The  peak  spreading  caused  by  elec- 
tronic noise  can  be  neglected  in  spectrum  analysis  if  the  same 
equipment  is  used  to  obtain  both  the  reference  peak  and  the  data 
to  be  analyzed. 

Relating  Observed  Counts  To  Activity 

The  pulse-height  spectra  obtained  from  alpha  particles 
emitted  by  isotopes  of  the  actinide  series  are  expected  to  be 
a composite  of  peaks  described  above.  In  order  to  analyze  these 
alpha  spectra,  it  is  necessary  to  determine  the  area  under  each 
peak  in  the  spectra.  After  the  area  has  been  determined,  it 
is  necessary  to  determine  the  activity  of  the  contributing  iso- 
topes. The  area,  which  is  the  total  number  of  counts  observed 
for  each  isotope,  can  be  related  to  the  activity  of  the  sample 
by  the  following  equation: 

Nj./t„  = AxGxf  xf  x 6 x f.  x f (6) 

tc  SSW  toc 


N.  = total  number  of  counts  = the  total 
area  for  each  isotope 

t = counting  time  (sec) 
c 

A = activity  of  source  in  nuclear 

transformations  per  second  (ntps) 

G = solid  angle  subtended  by  the  sensitive 
region  of  the  detector  relative  to  the 
source 

fflS  = scatter  of  radiation  from  surroundings 
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f = absorption  by  air,  detector  window, 
w and  dead  layer 

6 = intrinsic  efficiency  of  detector  = 

fraction  of  particles  which  enter 
sensitive  volume  and  produce  a pulse, 
assuming  no  dead  time  {#cnts/#a's) 

fj.  = correction  for  dead  time 

f = number  of  alpha  particles  observed  for 
a each  nuclear  transformation  (#a's/nt) 

If  it  is  assumed  that  equipment  is  not  changed,  scattering 
does  not  occur,  and  dead  time  does  not  exist,  Eq  (6)  can  be 
simplified  and  can  be  solved  for  the  activity  as  follows: 

A = (Nt/tc)  x (l/fa)  x (1/G)  x (1/6)  (7) 

Therefore,  if  the  area  of  each  peak  in  the  spectrum  can 
be  calculated,  then  the  activity  of  an  unknown  sample  giving 
rise  to  those  peaks  can  be  determined. 

Program  ALPHAFIT 

Program  AUTOFIT  was  written  to  analyze  charged  particle 
spectra  from  a mass  spectrograph.  The  modifications  on  AUTOFIT 
included  (1)  increasing  the  size  of  the  reference  peak  from 
I4.O  channels  to  2048  channels,  (2)  disposing  of  the  peak  search 

■ 

routine  and  Q-value  options,  and  (3)  decreasing  the  number  of 
peaks  to  be  analyzed  from  200  peaks  to  20  peaks.  The  peak  fit- 
ting algorithms  were  not  modified.  After  the  above  modifications 
were  performed,  the  name  of  the  program  was  changed  to  ALPHAFIT 
to  avoid  confusion.  A detailed  description  of  the  peak  analysis 
is  listed  below. 

13 
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Detailed  Peak  Analysis . A complex  spectrum  can  be  repre- 
sented by  an  array  [E]  which  is  formed  from  the  sum  of  a set 
of  resolved  components  Z(x).  Each  element  of  the  array  can  be 
expressed  as 

L 

E^  = ^2  ai^i (x^)  (6) 

where  Z^(x^)  designates  the  number  of  counts  in  the  ith  channel 
of  an  unnormalized  reference  peak  located  at  an  abscissa  x^,  and 
a^  is  a normalization  factor.  The  sum  extends  over  L individual 
components  in  the  complex  spectrum.  The  quantity  Zi(x^)  is  de- 
fined as  zero  if  the  reference  peak  located  at  x^  does  not  ex- 
tend to  the  ith  channel.  The  quantities  x^  and  a^  are  varied 
until  the  array  [E]  is  closely  matched  to  the  array  of  ordinates 
[Y]  of  the  data.  This  is  quantitatively  accomplished  by  forming 
a chi-squared  function 

N P P 

P = I (Y.-E. ) / ( 6Y. ) (9) 

i=1  1 1 1 

summed  over  all  the  channels  of  the  data.  The  function  F is 
then  minimized  with  respect  to  a^  and  x^.  During  minimization, 
the  weights  6Y^  are  taken  to  be  V Y^IO  , where  Y^  is  the  ordi- 
nate of  the  ith  data  point.  Since  Poisson  Statistics  apply  in 
counting  experiments,  this  weighting  factor  is  approximately 
equal  to  the  square  root  of  the  variance.  The  addition  of  10 
to  Y^  slightly  increases  the  relative  weighting  of  the  data  near 
the  tops  of  the  peaks  and  prevents  division  by  zero. 
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The  chi-squared  function  represents  the  ratio  of  the  ob- 
served spread  between  the  data  and  the  reference  peak  to  the 
expected  spread  between  the  data  and  the  reference  peak.  By 
minimizing  the  chi-squared  function,  the  best  approximation  to 
the  real  data  is  obtained. 

Program  ALPHAPIT  separates  the  minimization  of  F with 
respect  to  the  a^'s  and  the  x^'s.  If  the  x^'s  are  known  from 
the  initial  estimates  of  peak  positions,  the  a^'s  can  be  ob- 
tained by  setting  the  partial  derivative  of  F with  respect  to 
a1  equal  to  zero.  The  set  of  linear  equations,  which  need  to 
be  solved  in  order  to  find  the  correct  value  of  a^  and  which 
will  allow  minimization  of  F,  is  set  up  in  Subroutine  FCN. 

Standard  matrix-inversion  techniques  are  then  used  to  solve 
these  equations. 

The  position  of  the  peaks  are  established  by  a variable 
metric  minimization  procedure  (Ref  7:ANL-5990)  coded  in  a set 
of  subroutines.  As  the  x^'s  are  varied  by  Subroutine  SHIFT, 
the  function  F is  recalculated  (each  time  with  a best  set  of 
a^'s)  and  inspected.  If  the  change  in  F between  successive 
complete  iterations  is  less  than  a convergence  criterion,  the 
procedure  is  terminated.  A final  best  set  of  a^’s  is  calculated, 
and  the  results  are  transmitted  to  the  calling  program. 


«**•***, 
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III.  Experimental  Equipment  and  Procedures 
Introduction 

In  order  to  test  Program  ALPHAFIT,  it  was  necessary  to 
obtain  selected  alpha  spectra  with  certain  parameters  fixed. 

For  this  reason  experiments  were  performed  on  site  with  existing 
alpha  sources  such  as  Am2^1,  Po210,  and  Pu2^0. 

Equipment 

. 

The  alpha  spectroscopy  was  performed  using  a semiconductor 
detector  in  an  evacuated  chamber  to  maximize  the  resolution 
capabilities  of  the  system. 

Vacuum  Chamber.  The  vacuum  chamber  was  a two-piece,  cylin- 
drical steel  canister.  The  top  of  the  canister  was  fastened  to 
a wood  support  which  held  the  chamber  assembly  approximately 
one  foot  above  the  table  top.  Three  openings  were  located  in 
the  top  of  the  canister,  two  of  which  were  an  inlet  and  an  out- 
let for  the  vacuum  pump.  The  remaining  opening  allowed  the 
direct  connection  of  the  preamplifier  to  the  detector  by  two 
BNC  connectors.  The  detector  was  mounted  to  the  top  of  the  can- 
ister and  was  facing  downward  toward  the  moveable  source  plat- 
form. The  bottom  portion  of  the  canister  was  removeable  to 
allow  the  changing  of  sources  and  was  held  in  place  by  the  vacuum 
created  by  the  roughing  pump. 


Detectors.  A high  resolution  detector  with  an  active  area 
2 

of  25  mm  was  used  for  data  acquisition.  The  small  surface  of 
this  detector  required  a necessary  increase  in  the  counting 


times  and  required  the  reduction  of  the  source  to  detector 
distance.  The  resolution  obtained  by  this  detector  was  l8  keV 
FWHM  for  the  5*486  MeV  peak  of  Am2*4-1  with  a source  to  detector 
distance  of  1 cm.  The  resolution  was  not  appreciably  affected 
by  variations  in  source  to  detector  distances  between  1 to  3 cm; 
however,  the  counting  times  for  the  sources  used  were  much  too 
long  for  distances  above  1 cm. 

Amplifiers . An  ORTEC,  Model  142-B,  preamplifier  was  used 
in  conjunction  with  a linear  amplifier  and  a biased  amplifier 
to  obtain  the  desired  alpha  spectra.  The  linear  amplifier  was 
a Tennelec,  Model  TC203  BLR,  and  the  biased  amplifier  used  was 
a Tennelec,  Model  TC250A. 

Experimental 

Throughout  the  entire  set  of  experiments,  the  maximum 
amount  of  separation  between  the  peaks  of  the  spectra  was  de- 
sired. This  required  using  a biased  amplifier  which  would  cut 
out  signals  of  low  energy  and  expand  the  remaining  signals  to 
allow  the  maximum  peak  spread  possible. 

Energy  Calibration.  All  energy  calibrations  were  performed 
using  the  5*4-86  MeV  and  the  5*443  MeV  peaks  of  Am2^  . The  energy 
calibration  obtained  from  Am2^  was  used  when  sources,  such  as 
Po^1®  (which  ha3  only  one  peak)  and  Pu^°  (which  had  low  activity), 
were  used.  It  was  assumed  that  the  energy  calibration  obtained 
from  Am2^1  could  be  used  with  other  sources  as  long  as  the  linear 
amplifier  and  the  biased  amplifier  remained  at  the  same  settings. 
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Sources  and  Source  Preparation.  During  the  course  of 
this  study,  six  different  sources  were  used.  Three  of  the 
sources  had  been  prepared  at  an  earlier  date.  These  sources 
were  Am2^1  with  an  activity  of  0.1  (iCi  and  two  plutonium  sam- 
ples labeled  Pu2^-0  - 1 and  pu2^0  - 2 with  activities  of 
2.583  x 10*  cc's/min  ± 1.5%  and  2.598  x 10^  oc's/min  + 2%, 
respectively. 

Two  polonium  sources  were  made  during  the  course  of  this 

210 

study.  The  first  polonium  source,  labeled  Po  - 1,  was  made 

210 

by  the  electrochemical  displacement  of  Po  onto  a silver  disk 
from  a RaDEF  solution.  The  plating  occurred  because  of  the 
following  reaction: 

Po+++  + 3Ag  - Po  + 3Ag+ 

This  disposition  occurs  without  the  application  of  external 
current  when  the  RaDEF  is  dissolved  in  nitric  acid  (Ref  6:100). 
The  activity  of  this  polonium  source  was  roughly  approximated 
to  be  of  the  order  of  0.03  jiCi. 

The  second  polonium  source  was  prepared  using  a method 
described  by  H.  W.  Kirby  (Ref  3:2043-2047).  A half  inch  plat- 
inum square  was  rinsed  in  water;  then  fired  to  a red  heat  for 
ten  to  twenty  seconds.  After  cooling,  the  square  was  sprayed 
with  Krylon  with  the  center  of  the  square  masked  by  a brass  rod 
^ approximately  7 mm  in  diameter.  The  RaDEF  solution  which  was 

in  0.8n  HNO-j  was  dried  under  a heat  lamp;  then  redissolved  with 
1.5N  HCL.  The  solution  was  placed  in  the  center  of  the  platinum 


18 


square  which  was  then  placed  on  a hot  plate  until  the  solution 
was  completely  evaporated.  Then,  3N  acetic  acid  was  added  to 
the  center  of  the  square  and  allowed  to  dry.  The  source  was 
rinsed  with  water  and  then  acetone  to  remove  the  Krylon.  The 
activity  was  calculated  to  be  approximately  1 pCi. 


IV 


Data  Reduction 


Three  series  of  tests  were  performed  on  Program  ALPHAFIT 
to  determine  its  peak  fitting  capability. 

Test  Series  I involved  spectra  composed  of  two  Gaussian 
shaped  peaks  with  statistical  fluctuations  which  were  generated 
to  emulate  spectra  obtained  by  experimental  methods.  The  pur- 
pose of  this  series  of  tests  was  to  determine  the  error  in- 
volved in  determining  the  areas  under  the  peaks  when  parameters 
such  as  peak  amplitude,  peak  separation,  and  magnitudes  of 
pulse-height  intervals  were  varied. 

Test  Series  II  involved  actual  experimental  spectra.  The 
emphasis  in  this  series  of  tests  was  to  extract  a single  alpha 
peak  shape  from  a complex  alpha  spectra  by  a method  developed 
by  P.  De  Regge  (Ref  2:269-280). 

Test  Series  III  involved  experimental  spectra  and  emphasis 
was  again  placed  on  obtaining  a single  alpha  peak  shape.  How- 
ever, in  this  test  series,  the  single  alpha  peak  was  obtained 
210 

from  Po  , a single  alpha  emitter. 

Test  Series  1^ 

In  this  test  series,  a Gaussian  pulse-height  spectrum  with 
statistical  fluctuations  was  used  to  generate  emulated  pulse- 
height  spectra.  Studies  were  made  of  the  capabilities  of  the 
computer  program  to  resolve  two  peaks  whose  relative  positions 
and  amplitudes  were  varied  along  with  variation  in  the  magni- 
tude of  the  pulse-height  interval  used  to  generate  the  simulated 
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spectra.  The  error  involved  in  each  study  was  reported  in  terms 


of  relative  percent  error  in  peak  area  (RPE^). 

RPEa  = J (computed  area  - actual  area)/actual  area J x 100  (10) 

If  it  is  assumed  that  the  statistical  fluctuations  in  counts 
are  described  by  a Poisson  distribution,  then  the  standard  de- 
viation in  the  number  of  counts  is  equal  to  the  square  root  of 
the  average  total  count  or  the  square  root  of  the  area  under  the 
peak.  In  this  report  the  criterion  that  the  area  computed  should 
be  within  three  standard  deviations  of  the  actual  area  was  used 
in  order  for  the  results  of  the  analysis  to  be  valid.  The  sat- 
isfaction of  this  criterion  would  allow  the  user  to  state,  with 
99.7%  confidence,  that  the  limiting  factor  in  the  analysis  is 
not  the  program  but  the  counting  statistics.  This  criterion  can 
be  expressed  in  terms  of  relative  percent  error  in  peak  area  due 
to  statistical  fluctuations  (RPEy-£- ). 

RPEVT  = (3V  actual  area/actual  area)  x 100  (11) 

The  results  of  the  tests  described  above  are  plotted  in 
Figures  I|.  and  5*  In  Figure  k the  relative  percent  error  in  peak 
area  is  plotted  versus  the  separation  of  the  peaks  in  units  of 
full  width  at  half  maximum,  whereas  in  Figure  5 the  relative 
percent  error  in  peak  area  is  plotted  versus  the  magnitude  of 
the  pulse-height  interval  in  units  of  o/channel  where  a is  the 
standard  deviation  of  the  Gaussian  distribution.  In  each  figure 
the  error  resulting  in  the  determination  of  Peak  2 is  being 
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plotted.  Peak  2 (P2)  being  to  the  left  of  Peak  1 (Pi)  in  the 
pulse-height  spectrum.  The  two  heavy  dashed  lines  running  hor- 
izontally across  each  graph  signifies  the  relative  percent  error 
in  peak  area  for  one  and  three  standard  deviations.  It  is  assumed 
that  the  results  shown  in  Figures  4 and  5 are  applicable  to  other 
complex  alpha  spectra  obtained  by  experimental  methods.  The  user 
must  make  certain  that  the  spectrum  to  be  analyzed  has  a peak 
separation  and  pulse-height  interval  which  will  allow  the  areas 
of  the  peaks  to  be  determined  within  three  standard  deviations. 
This  will  not  guarantee  that  the  areas  of  the  peaks  in  the  ex- 
perimental spectrum  will  be  determined  within  three  standard 
deviations,  but  affords  the  best  possible  chance  of  it  occurring. 

The  program  appears  to  operate  more  efficiently  as  the  peak 
separation  is  decreased  and  as  the  pulse-height  interval  is  in- 
creased. This  unrealistic  result  which  occurs  for  peaks  of  un- 
equal amplitude  is  suspected  to  result  from  the  minimization 
algorithm  used  by  the  program  (see  page  15).  The  exact  cause 
was  not  investigated  in  this  report;  therefore,  the  results 
plotted  in  Figures  4 and  5 will  be  used. 


22 


< PI /P2=1 0 


P1/P2-1 


Figure  5*  Vs*  Width  of  Pulse-height  IntervRl (o) 


Test  Series  II 


After  Program  ALPHAFIT’s  peak  fitting  capability  was  de- 
termined using  peaks  of  known  areas,  analysis  of  actual  alpha 
spectra  began. 

The  first  alpha  spectra  were  obtained  from  Am2^  . The 
peaks  resulting  from  Am2**1  had  a separation  greater  than  1 FWHM. 
If  it  is  valid  to  apply  the  results  obtained  in  the  limit  of 
precision  test  on  the  ideal  peaks  to  actual  alpha  spectra,  the 
area  under  the  peaks  should  be  able  to  be  calculated  within 
three  standard  deviations  if  the  correct  reference  peak  is  used. 
The  problem  which  resulted  was  how  to  obtain  the  proper  reference 
peak.  All  available  sources  of  sufficient  activities  emitted  at 
least  three  different  energy  alpha  particles.  Each  alpha  par- 
ticle creates  a peak  with  a long  exponential  tail.  When  more 
than  one  peak  exists  in  the  spectrum,  confusion  exists  as  to 
exactly  how  much  each  peak  contributes  to  the  tail. 

Since  it  was  impossible  to  obtain  a reference  peak  directly 
from  the  Am2**1,  a method  developed  by  P.  De  Regge  was  applied 
(Ref  2:269-280).  De  Regge  based  his  empirical  relation  on  ob- 
servations made  on  a single  alpha  peak.  It  was  suspected  that 

21 0 

the  single  alpha  peak  was  obtained  from  Po  . The  empirical 
relation  was  derived  from  an  alpha  peak  plotted  in  2048  channels 
on  semilog  paper  (schematically  shown  in  Figure  6).  The  low 
energy  side  of  the  peak  was  divided  into  three  regions.  The 
first  region  ranges  from  0 to  3 MeV  below  the  peak  energy  and 
can  be  represented  by  a constant  flat  line  (a-*). 
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Figure  6.  Schematic  of  De  Regge's  Alpha  Peak 

The  second  region  can  be  best  represented  by  a straight  line  (b-c) 
with  a slope  S.  The  third  region  can  best  be  represented  by  a 
straight  line  (d-e)  with  a slope  R,  which  is  much  greater  than 
slope  S.  Considering  the  spectrum  to  the  left  of  the  peak  value, 
the  following  equation  was  proposed: 

Ci+1  = C.eP(i)  02) 

where 

= content  of  channel  i 

F(i)  = a function  of  the  channel 
number  i 

S for  the  region  (b-c) 

R for  the  region  (d-e) 


The  transition  region  (c-d)  between  the  two  straight  line  sec- 
tions can  be  represented  by 


C<1)  - 


where  the  function  F(i)  is  best  fit  by  a hyperbolic  tangent 


function  of  the  form 


F(i)  = (R+S/2)  + [(R-S/2)tanha(i-k) J 


where  a is  a coefficient  to  accomodate  various  transitions  from 


the  tail  to  the  peak  regions.  Negative  as  well  as  positive 


values  are  required  for  the  argument  of  the  hyperbolic  tangent; 


therefore,  the  origin  of  the  i-axis  is  shifted  by  an  amount  k 


to  coincide  with  the  bending  point  of  the  hyperbolic  tangent. 


Integration  of  Eq(ll|)  and  substitution  into  Eq(13)  yields 


C(i)  = 21/Z  K[ezR(i-k)  + ezS(i-k)]l/z 


with  K being  the  constant  of  integration  and 


z m 2 o/(R-S) 


The  tail  region  (b-c)  is  characterized  by  channel  numbers 


i<xk.  The  channel  contents  in  this  region  yield  a straight  line 


with  slope  S on  a semilog  plot. 


lnC(i)  = InK  - l/zln2  + S(i-k) 


The  straight  line  with  slope  R,  fitting  the  logarithm  of 


the  channel  contents  in  the  region  (d-e),  is  given  by 


d 


lnC(i)  = InK  - l/zln2  + R(i-k) 


(17) 


The  intersection  point  of  these  straight  lines  is  located  at 
channel  number  k. 

Eq(l5)  is  fit  to  the  data  by  obtaining  CQ,  the  tail  height 
between  a and  b,  and  determining  the  slope  R from  a least  squares 
fit  of  Eq(l7)  to  the  logarithm  of  the  contents  of  the  channels 
between  d and  e.  The  slope  S was  found  to  be  constant  in  sev- 
eral spectra  and  had  a value  of  O.OOI4..  The  point  b where  the 
tail  region  changes  from  a constant  value  to  a gently  increasing 

curve  was  chosen  as  channel  number  (k-700)  according  to  obser- 

. 

vations  made  on  several  alpha  spectra.  Substitution  of  the  above 
into  Eq(l6)  yields 

lnCQ  = InK  - l/zln2  - 2.8  (18) 

which,  when  substituted  into  Eq(l7)»  yields 

lnC(i)  = lnCQ  + 2.8  + R(i-k)  (19) 

which  allows  the  calculation  of  the  location  of  point  k on  the 
abscissa.  The  value  of  z is  determined  only  by  the  ratio  a/R 
since  the  value  of  S is  much  smaller  than  R.  The  ratio  a/R  is 
expected  to  be  independent  of  the  resolution  and  of  the  shape  of 
the  alpha  peaks,  because  both  a and  R increase  with  better  reso- 
lution. Values  for  2a/R  have  been  found  to  be  nearly  constant 
and  equal  to  0.59  for  typical  resolutions  between  20  and  80  keV. 

The  value  of  K is  calculated  by  replacing  C(i)  in  Eq(l5) 
with  the  real  contents  of  a few  undisturbed  channels  near  the 
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top  of  the  peak  and  by  solving  the  equation  for  K.  The  mean 
value  of  K is  then  substituted  back  into  Eq(l5)  to  calculate 
the  channel  contents  in  the  region  (b-e).  The  high  energy 
edge  of  the  peak  is  assumed  undisturbed  and  is  simply  inte- 
grated. 

An  attempt  was  made  to  apply  the  method  described  above 

to  obtain  a reference  peak  from  an  Am^*^  spectrum.  This  was 

done  to  determine  whether  the  method  developed  by  De  Regge  was 
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applicable  to  alpha  spectra  other  than  Po  . However,  unlike 
De  Regge,  only  512  channels  were  used  to  store  the  spectrum. 
This  was  done  to  allow  for  easier  manipulation  of  data.  The 
Am^^  spectrum  was  obtained  in  5l2  channels  with  the  aid  of  a 
biased  amplifier  which  eliminated  the  low  energy  end  of  the 
spectrum  and  expanded  the  remaining  spectrum.  This  allowed  the 
prominent  alpha  peaks  to  be  well  separated  but  at  the  cost  of 
losing  the  constant  part  of  the  tail  (Region  (a-b),  Figure  6). 
It  was  impossible  to  determine  the  value  of  k using  the  method 
De  Regge  described  above;  therefore,  a simpler  method  was  de- 
vised. 

An  assumption  was  made  from  the  outset  that  the  addition 
of  a biased  amplifier  would  not  change  the  values  of  the  con- 
stants obtained  by  De  Regge.  A simplified  method  was  then 
developed  using  this  assumption  and  using  the  fact  that  the 
f - point  k is  located  at  the  intersection  of  the  lines  drawn  from 


(b— c)  and  (e-d)  (see  Figure  6).  The  line  (e-d)  and  slope  R 
were  easily  obtained  by  applying  the  method  of  least  squares 
to  the  contents  of  a few  channels  near  the  top  of  the  5 .14-86  MeV 


spectrum.  However,  line  (b-c)  was  not  so 


easily  obtained,  and  a few  more  assumptions  had  to  be  made. 

The  assumptions  made  in  obtaining  line  (b-c)  were 

a)  each  tail  could  be  approximated  by  an 
exponential, 

b)  the  slope  of  each  tail  is  the  same, 

c)  the  total  tail  is  composed  of  a linear 
combination  of  the  contribution  from 
each  peak, 

d)  each  peak  contributed  an  amount  to  the 
tail  proportional  to  the  area  under  its 
peak,  and 

e)  the  tail  contributions  above  the  S .i|-Q6>  MeV 
peak  are  negligible. 

The  total  tail  resulting  from  the  three  main  alpha  peaks  can 
be  represented  by 


Ni  = 


V 


-S(xp1-Xi) 


-s<vv 


C^e 


”S(Xp3”xi) 


(20) 


with 

= counts  in  the  total  tail  in  channel  i 

S = slope  of  tail  as  was  defined  by  De  Regge 

Xd1*  xd2*  xd3  = ciiannel  number  of  peak  1,  peak  2, 

P P Pi  peak  3,  respectively 

, C2»  Co  = relative  abundance  of  peak  1,  peak  2, 
peak  3,  respectively 

x^  a channel  number  of  interest 

Peak  1 = 5*486  MeV  peak  of  Am*^ 

Peak  2 * 5*443  MeV  peak  of  Am*^ 

Peak  3 - 5*389  MeV  peak  of  Am^1 


rr 

' 

The  slope  S was  calculated  by  using  a least  squares  fit 
on  the  tail  of  the  Am2^1  spectrum.  The  peak  locations  were 
known  and  the  constants  C1  and  C 2 could  be  written  in  terms  of 
C^.  The  constant  could  be  calculated  after  substituting  a 
known  value  for  N^,  along  with  the  corresponding  channel  num- 
ber, into  Eq(20).  With  known,  the  tail  resulting  only  from 
the  5*486  MeV  peak  could  be  calculated.  The  point  k was  visu- 
ally obtained  by  extending  the  line  formed  by  this  tail  until 
it  intersected  line  (d-e). 

Having  obtained  the  slopes  S and  R and  the  value  of  the 
point  k,  the  next  step  was  to  determine  the  value  of  the  inte- 
gration constant  K in  Eq(l5)*  The  value  of  K was  calculated 
by  substituting  a few  values  of  C(i)  in  Eq(l5)  with  the  real 
content  of  a few  undisturbed  channels  near  the  top  of  the 
5.486  MeV  peak.  The  mean  value  of  K was  substituted  back  into 
Eq(l5)  60  determine  the  intermediate  channel  contents  between 
the  tail  and  the  top  of  the  peak. 

The  values  resulting  from  Eq(l5)  were  checked  against  one 
important  criterion.  Since  the  actual  data  was  being  used  to 
determine  a reference  peak  shape,  it  would  be  necessary  to  re- 
quire that  the  contents  in  every  channel  calculated  by  Eq(l5) 
be  less  than  or  equal  to  the  contents  in  the  corresponding 
channel  of  the  real  spectrum.  If  the  calculated  peak  violated 
this  criterion,  the  peak  would  be  considered  invalid.  Several 
different  Am2^1  alpha  spectra  were  analyzed  using  a reference 
peak  obtained  from  Eq(l5)*  The  results  of  one  such  analysis 
are  listed  in  Table  I.  The  results  from  this  analysis  are  not 
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directly  comparable  to  those  obtained  by  De  Regge  because  his 
results  represent  an  average  of  13  separate  measurements.  Only 
a few  of  the  many  analyzes  performed  using  the  simplified  method 
were  able  to  satisfy  the  above  criterion  and  so  this  method  was 
not  used. 

Table  I 

Am^^  Analyzed  Using  Reference  Peak  Obtained  by 


Eq ( 1 5 ) and 

Simplified  De 

Regge  Method 

Alpha  I 

Alpha  II 

Alpha  III 

Peak  energy 

5.486  MeV 

5.443  MeV 

5.389  MeV 

Peak  area 

240799 

37749 

5187 

Branching  factor 

86.0% 

12.7% 

1.3% 

Computed  branching 
factor 

84.4% 

13.2% 

1.8% 

De  Regge' s computed 
branching  factor 

81 .8% 

14.5% 

3.7% 

♦ (From  Ref  4*435 ) 


Test  Series  III 

This  series  of  tests  involved  using  the  ability  of  Program 
ALPHAFIT  to  accept  a3  reference  peak  as  a table  of  values  spec- 
ifying the  peak  shape  channel  by  channel.  The  tests  series 

21 0 

evolved  around  the  single  alpha  emitter,  Po  . There  were 

21 0 

several  desired  objectives  in  using  Po  . These  objectives 
included  reproducing  De  Regge' s work,  with  hopes  of  obtaining 
new  constants,  which  would  allow  more  reliable  reproduction  of 
the  reference  peak  in  different  spectra  and  evaluating  the  change 

31 


of  the  resulting  peak  shape  with  source  thickness.  However, 

before  any  of  these  steps  were  begun,  a series  of  tests  were 

210 

performed  in  which  an  actual  Po  spectrum  was  used  as  a 

reference  peak.  The  empirical  expression,  developed  by  De 

Regge,  was  an  approximation  to  a single  alpha  peak  such  as 
2l0 

obtained  from  Po  . The  best  fit  that  could  be  obtained 

would  result  when  the  mathematical  expression  exactly  described 

the  single  alpha  peak.  De  Regge  states  that  alpha  peaks  have 

the  same  shape  independent  of  their  energy  when  thin  sources 

are  used  (Ref  2:271).  If  this  assumption  is  viable,  then  the 

limit  of  precision  of  the  procedure  developed  by  De  Regge  can 

be  quickly  determined  by  using  the  single  alpha  peak  (such  as 

obtained  from  Po  ) as  the  reference  peak  in  the  analysis. 

The  results  of  this  analysis  are  listed  in  Table  II.  Using 
210 

Po  directly  as  a reference  peak  did  not  seen  to  give  better 

results  than  those  obtained  by  De  Regge' s method.  The  larger 

error  in  branching  ratios  could  have  resulted  from  the  effects 

of  varying  source  thicknesses.  The  resulting  error  could  mesn 

that  the  sources  are  of  different  thicknesses  and  do  not  qualify 

as  a thin  source;  hence,  the  sources  are  not  expected  to  be  fit 

2 io 

using  Po  . The  variations  in  source  thicknesses  are  clearly 
evident  in  the  analysis  of  Po210  - 2 using  Po210  - 1 as  a ref- 
erence peak.  The  residual  is  very  large  and  is  due  to  a dif- 
ferent source  thickness  and  a different  level  of  source  contam- 
^ ination.  A much  improved  analysis  is  expected  to  result  when 

the  sources  are  nearly  the  same  thickness. 
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Table  II 

2 10 

Isotopes  Analyzed  Using  Po  - 1 Directly 
as  the  Reference  Peak 

— 

Po210  - 1 Alpha  I 

Energy  (MeV)  5»305 

Computed  Peak  Area  = 302894 


Actual  Peak  Area  =■  302946 

Residual  = 5 2 RPEA  = 0.02% 

3lfT  = 1651  RPE  ^ = 0.55% 


to2*1 

Alpha  I 

Alpha  II 

Alpha  III 

Energy  (MeV) 

5.436 

5.443 

5.389 

Computed  Peak  Area 

811711 

80328 

1178 

Branching  Factor 

86% 

12.7% 

1.3% 

Computed  Branching 
Factor 

90.9% 

8.9% 

0.13% 

Computed  Total  Area 

= 893217 

Actual  Total  Area 

= 896127 

Residual 

= -4910 

Rpea  = 0.55% 

3 V"a 

= 2843 

RPE^  = 0.32% 

T&°  - 1 

Alpha  I 

Alpha  II 

Alpha  III  Alpha  IV 

Energy  (MeV) 

5.499 

5-456 

5.168  5.123 

Computed  Peak  Area 

21939 

8928 

702697  493291 

Branching  Factor 

0.1224% 

0.0476% 

74.75%  23.6% 

Computed  Branching 
Factor 

1.8% 

0.73% 

57.3%  40.2% 
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Table  II  (cont.) 


Pu2l+0  - 1 (cont.) 


Computed  Total  Area 
Actual  Total  Area 
Residual 


1229505 

1226633 


2872 

3323 


rpea  = 0.23% 

RPE,/r  = 0.27% 


Pu2^0  - 2 

Alpha  I 

Alpha  II 

Alpha  III  Alpha  IV 

Energy  (MeV) 

5-499 

5-456 

5.168 

5.123 

Computed  Peak  Area 

2501 

514 

90664 

30760 

Branching  Factor 

0.1224% 

0.0476% 

74.75% 

23.6% 

Computed  Branching 
Factor 

2% 

0.41% 

72.9% 

24.7% 

Computed  Total  Area 

= 127271 

Actual  Total  Area 

= 125341 

Residual 

= 1930 

rpea  = 

1.54% 

3\rA 

= 1062 

RPE,s  * 

0.85% 

Po210  - 2 

Alpha  I 

Energy  (MeV) 

5.305 

Computed  Peak  Area  = 

8745688 

Actual  Peak  Area  - 

8761782 

Residual  = 

- 1 6094 

RPEa  = 

0.18% 

3 V~a 

8880 

rpeVa  ■ 

0.10% 
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Table  II  (cont.) 


Composite  Spectrum 

(PO210  - 2, 

Am 21+1 , 

Pu 

24°  _ 2) 

Am21*’ 

Alpha  I 

Alpha 

II 

Alpha  III 

Energy  (MeV) 

5.486 

5.433 

5.389 

Computed  Peak  Area 

853899 

80306 

946 

Branching  Factor 

86% 

12.7% 

1.3% 

Computed  Branching 
Factor 

91.3% 

8.6% 

0.1% 

Computed  Total  Area 

= 935151 

Actual  Total  Area 

= 898^27 

Residual 

= 37024 

rpea  = 4.12% 

3Va 

= 2843 

RPEVa  = °*32* 

Composite  Spectrum 

(cont. ) 

Pu21*0  - 2 

Alpha  I 

Alpha 

II 

Alpha  III  Alpha  IV 

Energy  (MeV) 

5.499 

5.456 

5.168  5.123 

Computed  Peak  Area 

0 

0 

60259  15170 

Branching  Factor 

0.1224% 

0.0476% 

74.75%  23.6% 

Computed  Branching 
Factor 

0% 

0% 

79.9%  20.1% 

Computed  Total  Area 

= 75429 

Actual  Total  Area 

= 125341 

Residual 

= -49912 

rpea  = 39.82% 

3\rA 

= 1062 

r?ev£  = 0.85% 
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Table  II  (cont.) 


Composite  Spectrum 

(cont. ) 

Po210  - 2 

Alpha  I 

Energy  (MeV) 

5.305 

Computed  Peak  Area 

* 1499601*. 

Actual  Peak  Area 

= 1810354 

Residual 

* -310750 

RPEa  = 17.17% 

3V"A 

= 4036 

RPE  = 1.30% 
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The  last  set  of  tests  in  Test  Series  III  involved  using 
an  entire  spectrum  of  an  isotope  as  a reference  peak.  The 
composite  spectrum  to  be  analyzed  is  known  to  contain  certain 
isotopes  (Pu2^0  - 2,  Po210  - 2,  Am2*4-1 ) . Several  experiments 
were  performed  in  which  the  spectrum  from  each  isotope  was 
obtained.  The  three  spectra  were  then  added  together  to  simu- 
late a composite  spectrum  containing  the  three  isotopes.  During 
the  analysis,  the  spectrum  from  each  isotope  was  used  as  the 
reference  peak  for  that  isotope.  The  results  of  this  analysis 
are  listed  in  Table  III. 


Table  III 


Composite  Spectrum  Analyzed  Using  Respective  Spectrum 
from  each  Isotope  as  the  Reference  Peak 


Composite  Spectrum 

(Po210  - 2 

, Am2i;1, 

Pu21*0  - 2) 

Computed 

Actual 

Area 

Area 

Residual 

3\Ta 

Am21*-1 

2712695 

898127 

1814568 

2843 

Pu2^0  - 2 

2578378 

125341 

2453037 

1062 

Po210  - 2 

1995679 

1810354 

185325 

4036 

The  residuals  resulting  from  this  analysis  are  well  above  the  de- 
sired limit.  This  error  is  suspected  to  have  resulted  from  not 
utilizing  the  program's  full  capability.  The  analysis  was  per- 
formed as  follows: 
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a)  Spectrum  to  be  analyzed  read  into  program 

b)  Reference  peak  (entire  spectrum  of  one 
isotope)  read  into  program 

c)  Reference  peak  fit  to  the  data  and  results 
printed  out 

d)  Program  cleared 

e)  Spectrum  to  be  analyzed  read  back  into 
program 

f)  Next  reference  peak  read  into  program 

g)  Reference  peak  fit  to  the  data  and  results 
printed  out 

h)  etc. 


A better  analysis  would  result  if  after  the  spectrum  to 
be  analyzed  has  been  read  into  the  program  all  three  reference 
peaks  were  read  into  the  program  simultaneously.  The  program, 
working  up  to  its  full  potential,  could  then  vary  the  relative 
positions  and  amplitudes  of  the  reference  peaks  to  obtain  the 
best  fit  possible. 
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r.  Conclusions  and  Recommendations 

Throughout  this  thesis,  reference  was  made  to  the  "proper 
reference  peak".  The  "proper  reference  peak"  is  one  which  best 
represents  the  parent  population.  The  "proper  reference  peak" 
was  used  in  Test  Series  I.  The  spectra  in  these  tests  were  the 
parent  populations,  statistically  varied  to  simulate  actual  ex- 
perimental spectra.  However,  this  was  done  with  the  qualifica- 
tion that  both  the  reference  peak  and  the  spectra  were  from  the 
same  type  of  source,  having  the  same  thickness  and  backing  ma- 
terial. 

The  ability  of  Program  ALPHAPIT  to  fit  unknown  alpha  spec- 
tra depends  upon  the  ability  of  the  user  to  find  the  "proper 
reference  peak".  In  order  to  find  the  "proper  reference  peak", 

one  needs  to  study  a single  undisturbed  alpha  peak.  This  was 

2 1 0 

done  in  the  past  using  Po  (Ref  2:269-280).  Even  after  using 
the  single  alpha  peak,  De  Regge  had  to  try  to  find  some  empir- 
ical function  which  would  allow  him  to  approximate  part  of  the 
peak.  This  approach  gives  good  results,  but  errors  result  due 
to  the  assumptions  made  by  De  Regge.  Most  of  the  error  arose 
from  the  empirical  function  used  to  represent  the  peak.  Some 
error  could  have  arisen  from  the  assumption  that  a complex 
alpha  spectrum  is  composed  of  individual  alpha  peaks  which  do 
not  interfere  with  each  other,  but  only  add  linearly  to  one 
another.  With  the  aid  of  Program  ALPHAPIT,  the  approximations 
used  in  trying  to  find  an  analytical  function  to  describe  the 
alpha  peak  are  circumvented.  Instead,  the  single  alpha  peak 
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from  Po  can  be  directly  used  by  the  program  to  analyze  com- 
plex alpha  spectrum.  However,  the  assumption  that  complex 


spectra  can  be  linearly  composed  of  single  alpha  peaks  must 

still  be  used.  If  this  is  true  and  if  it  can  be  assumed  for 

the  moment  that  for  sufficiently  thin  sources  (Ref  2:269-200) 

the  alpha  peaks  will  look  the  same,  then  ALPHAPIT  should  give 
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the  most  precise  results  obtainable  assuming  the  Po  alpha 
peak  is  obtained  by  the  same  detector. 

The  last  series  of  tests  performed  were  bas*u  on  the  use 
of  an  entire  spectrum  as  a reference  peak.  This  would  allow 
circumvention  of  both  approximations  made  by  De  Regge.  No 
empirical  function  would  be  needed,  nor  would  there  be  a need 
for  the  assumption  of  linear  addition  of  single  alpha  peaks. 
This  should  allow  the  most  precise  analysis  possible.  The  fact 
that  large  errors  resulted  in  the  first  attempt  to  apply  this 
method  should  not  hinder  its  further  development.  Perhaps  a 
spectrum  stripping  approach  would  be  advisable  in  this  case. 

The  application  of  either  method  requires  known  samples 
of  high  activity  to  obtain  good  statistics  in  the  tail  region. 
If  the  same  isotope  is  to  be  used  to  obtain  this  high  activity 
source,  it  will  require  either  more  area  or  a thicker  source. 
This  is  expected  to  reduce  the  effectiveness  of  the  reference 
peak,  because  the  peak  shape  will  change.  The  change  in  peak 
shape  should  be  relatively  unchanged  if  thin  sources  are  used 
(Ref  3:20i|.3-2047)  • The  exact  error  incurred  would  require 
further  tests  using  varying  source  thicknesses. 

At  this  time,  there  appear  to  be  three  different  methods 
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which  would  allow  a possible  solution  to  the  alpha  fitting 
problem: 


I.  Reproduce  De  Regge ' s work  to  determine  the 
change  in  peak  shape  with  source  thickness. 

II.  Investigate  further  the  possibilityQof  using 
a single  alpha  emitter  such  as  Po^  as  a 
reference  peak. 

III.  Investigate  further  the  possibility  of  using 
each  respective  spectrum  in  the  composite 
spectrum  as  a reference  peak. 

It  would  be  advantageous  to  include  a correction  factor 
for  source  thickness  in  each  of  these  methods.  This  could  be 
accomplished  with  a library  of  reference  peaks  for  sources  of 
various  thicknesses.  An  example  of  such  a library  could  be 
taken  from  data  such  as  is  found  in  Figure  2.  Knowing  a few 
parameters,  such  as  the  shape  of  the  high  energy  edge  of  the 
peak  or  the  full  width  at  half  maximum,  would  allow  the  user 
to  enter  the  library  and  pick  the  appropriate  reference  peak. 
If  a sufficient  number  of  reference  peaks  are  catalogued,  lin- 
ear interpolation  could  be  assumed  between  reference  peaks. 
This  would  provide  an  unlimited  number  of  possible  reference 
peaks  for  a variety  of  source  thicknesses. 
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APPENDIX  A 


Program  ALPHAFIT 

Main  Program.  The  main  program  controls  the  reading  of 
the  input  cards,  initializes  variables,  checks  for  error 
conditions,  and  prepares  the  data  for  analysis. 

All  of  the  data  is  read  in  the  main  program  except  for 
the  spectrum  data  set,  which  is  obtained  from  Subroutine  SPCTRM. 
A subroutine,  named  ARYODR,  is  called  to  rearrange  background 
points  and  peak  positions  into  ascending  order  of  abscissas. 

As  the  data  is  prepared  for  analysis,  the  main  program 
establishes  the  logic  for  selecting  a reference  peak,  estimating 
peak  positions,  and  removing  background.  Two  subroutines,  named 
BGRND  and  REFPEK,  are  called  for  establishing  the  background  and 
the  reference  peak  shape,  respectively.  Afterwards,  a new  spec- 
trum with  background  subtracted  is  computed  with  the  requirement 
that  none  of  the  computed  ordinates  be  less  than  zero. 

Before  analysis  of  the  peaks  begins,  a final  check  is  made 
to  determine  whether  the  estimated  peak  locations  are  in  the 
range  of  the  spectrum.  Peaks  outside  the  range  are  deleted  from 
further  processing. 

Subroutine  RESLTS.  Subroutine  RESLTS  is  then  called  to 
process  the  data.  Subroutine  RESLTS  transmits  the  spectrum  to 
a subroutine,  called  PREPAR,  which  establishes  the  initial  error 
matrix,  the  weights  for  the  data  points,  and  the  convergence 
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criterion.  On  returning  from  Subroutine  PREPAR,  the  normalized 
fitted  peaks  are  summed  to  obtain  the  areas.  A printout  of  the 
channel  number,  the  raw  data,  the  background,  the  normalized 
individual  peaks  (with  background  included),  and  the  composite 
spectrum  of  the  fitted  peaks  can  be  obtained  by  the  user  upon 
request. 

A final  summary  page  lists  the  information  characterizing 
the  spectrum  and  tabulates: 

a)  peak  position  (channel  number), 

b)  error  in  peak  position  (channels), 

c)  number  of  counts  in  the  peak, 

d)  statistical  error  in  the  number  of  counts 
computed  by  the  program  during  analysis 
procedures, 

e)  sum  of  the  background  underneath  the  peak, 
and 

f)  estimated  error  in  the  number  of  counts 
(this  is  the  square  root  of  the  sum  of 
the  counts  and  the  background). 

Subroutine  REFPEK.  The  one  characteristic  that  made  this 
program  look  promising  was  the  ability  to  input  any  peak  shape. 
This  reference  peak  is  not  obtained  from  an  analytical  function 
but  is  expressed  as  a table  of  values  specifying  the  peak  shape, 
channel  by  channel,  with  linear  interpolation  assumed  between 
channels . 

The  reference  peak  is  read  into  the  computer  in  the  main 


program,  and  then  evaluated  by  Subroutine  REFPEK.  This  sub- 
routine uses  parameters,  such  as  peak  width  and  shape,  to 


Subroutine  BGRND.  Three  options  are  available  for  re- 
moving backgrounds  from  the  spectrum.  One  option  is  completely 
automatic;  the  other  two  require  information  from  the  user. 

Background  option  I allows  the  user  to  specify  coordinates 
of  a set  of  background  points  throughout  the  spectrum.  Sub- 
routine BGRND  then  calculates  a continuous  line  that  passes 
through  all  of  these  points.  The  regions  between  the  points 
are  calculated  using  Lagrange's  three-point  interpolation  for- 
mula. Three  successive  points  are  used  to  calculate  a curve 
between  the  first  two  points.  The  first  point  is  then  dropped, 
and  the  next  three  successive  points  are  used  to  calculate  the 
next  section.  This  results  in  a final  curve  which  is  a patch- 
work  of  parabolas.  The  computed  background  line  is  continuous, 
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but  its  derivative  may  be  discontinuous  at  each  background 
point.  The  background  is  never  allowed  to  be  less  than  zero. 

With  background  option  II,  the  spectrum  can  be  divided 
into  sections  with  a constant  background  value  in  each  section. 
The  limits  of  the  sections  and  the  background  value  of  those 
sections  are  supplied  by  the  user.  This  option  is  useful  when 
the  user  desires  to  have  a constant  background  for  the  entire 
spectrum. 

Background  option  III,  the  automatic  background  option, 
requires  that  the  estimated  locations  of  the  peaks  be  supplied. 
Subroutine  BGRND  examines  the  regions  between  such  peaks  and 
excludes  the  channels  that  contain  the  tails  of  the  peaks.  If 
such  a region  is  longer  than  ten  channels,  it  is  divided  into 
subregions  of  about  ten  channels  each.  The  abscissa  of  a back- 
ground point  is  taken  to  be  the  midpoint  of  the  abscissas  of 
the  subregion,  and  the  ordinate  of  a background  point  is  taken 
to  be  the  average  of  the  ordinates  of  the  data  points  in  the 
subregion.  Since  most  physical  spectra  have  backgrounds  that 
decrease  as  the  channel  number  increases,  a constraint  is  im- 
posed on  the  background  points.  The  ordinate  of  each  new  back- 
ground point  Y,  . is  compared  with  the  ordinate  of  the  last 

D,  1 

background  point  If  the  new  background  point  exceeds 

the  value  i-l+  2\  ^,i-l  * the  P°int  is  rejected.  After  at 
most  fifty  background  points  have  been  selected,  a subroutine, 
named  LSQPOL,  applies  the  method  of  least  squares  to  compute  a 
polynomial  of  order  less  than  four  through  the  background  points. 
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APPENDIX  B 


INPUT  CARDS  FOR  PROGRAK  ALPHAFIT 


Column  Variable 


Value 


Usage 


CARD  1 . Title  Card 

TITLE  (1),  I = 1,  8 
FORMAT  (8A10) 

TITLE  = Any  alphanumeric  string. 


1-80 

NOTE: 


If  columns  1-4  are  blank,  the  program  terminates, 
program  al3o  terminates  on  an  end-of-file  card. 


The 


CARD  2a.  Option  Control  Card 

NBACK,  NPKRD,  NREF,  NSSW1 , NSSW2,  MAXREF 
FORMAT  (14,  412,  14) 

NBACK  = 0 


1-4 


NOTE  2: 


5-6 


+n 


-n 


NOTE  1 : n < 50 . 


Background  determined  auto- 
matically by  program. 

Background  option  I.  The  X 
and  Y coordinates  for  n back- 
ground points  will  be  read 
in;  program  calculates  a 
curve  passing  through  all 
points. 

Background  option  II.  n 
background  intervals  will 
be  read  in  (background  is 
a constant  between  the  limits 
of  each  interval). 

If  n > 50,  the  program  will 
truncate  it  to  50. 


n > 3 for  background  option  I (positive  sign).  But 
also  see  Note  2 for  card  3« 


NPKRD 


+1 


Initial  peak  positions  to  be 
read,  as  values  of  channel 
number. 
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7-8 


NREF 


10 

12 

13-16 

CARD  2b. 

1 

5 

CARD  3. 

1-5 


NSSW1 


NSSW2 


MAXREF 


0 Reference  peak  to  be  se- 

lected from  spectrum. 

(This  option  is  generally 
not  recommended.) 

-1  Reference  peak  values  to 

be  read  from  separate  cards. 

0 Suppress  printed  output 

from  variable-metric-mini- 
mization routines. 

(Standard ( ) 

1 Print  the  v.m.m.  calculations. 

0 Suppress  printed  output  of 

raw  data,  resolved  peaks, 
and  composite  spectra. 

1 Print  the  above. 

The  maximum  number  of  channels 
in  the  reference  peak. 


Option  Control  Card 
KEPREF,  KEPDAT 
FORMAT  (II,  3 x,  II) 
KEPREF  * 0 

= 1 

KEPDAT  = 0 

= 1 


New  reference  peak  used  for 
each  spectrum  analyzed. 

Allows  the  reference  peak 
to  be  kept  for  next  analysis. 

New- data  spectrum  to  be 
analyzed. 

Allows  the  data  spectrum  to 
be  kept  for  next  analysis. 


Background  Data  Cards  (Not  used  for  NBACK  = 0.) 

DB ( I ) , BACK(I),  1=1,  NBACK 
FORMAT  (16F5.0) 

DB ( I ) = For  NBACK  > 0 (option  I): 
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11-l5»  etc 


6-10 


NOTE  1 

NOTE  2 


CARD  k 


1 


I 


DB(I)  is  the  channel 
number  at  which  a background 
value  is  to  be  specified. 

For  NBACK  < 0 (option  II): 

DB(I)  is  the  channel 
number  at  the  lower  end  of 
the  Ith  interval  for  which 
a background  value  is  to  be 
specified.  DB(1)  is  supplied 
by  the  program  and  may  be 
left  blank. 

BACK(I)  = The  value  of  the  background 

for  the  Ith  background  posi- 
tion (NBACK  >0),  or  for  the 
Ith  background  interval 
INBACK  <0 ) . 

The  pairs  DB(I),  BACK(I)  may  be  given  in  any  order. 

The  program  will  rearrange  them  in  increasing  order 
of  DB ( I ) . 

For  NBACK  > 0,  the  program  requires  backgrounds  for 
the  first  and  last  channels  of  data.  If  not  given 
explicitly,  the  program  will  use  BACK  (1)  for  the 
first  channel,  and  BACK  (NBACK)  for  the  last  channel. 
These  program-supplied  values  may  be  counted  for 
meeting  the  minimum  requirement  of  3 background  points. 


. Reference  Peak  Values  Card  (Use  only  for  NREF  < -1.) 

N 

FORMAT  (14) 

N * The  number  of  channels  in 

the  reference  peak. 

DUMCHN,  (ST0REF( J+K) , K = 1,  10) 

FORMAT  (F6.0/10(F6.0,  2x)) 

DUMCHN  = Every  tenth  channel  number 

of  the  reference  peak.  These 
values  are  not  needed  in  the 
reference  peak  and  so  are  not 
stored.  This  allows  cards 
punched  from  MCA  to  be  read 
without  conversion. 


k9 


t 


4 


1-6 
9-14 
etc . 


STOREF( J+K) 


The  number  of  counts  in  the 
( J+K ) th  channel  of  the  ref- 
erence peak,  reading  from 
the  back  edge  to  the  front 
edge. 


CARD  4b.  Reference  Peak  Parameter  Card  (Use  only  if  NREP  = 0.) 
STOIPR,  DELTA 1 , DELTA2,  DELTA3 
FORMAT  (4F10.5) 


STOIPR 


DELTA1 


DELTA2 


DELTA3 


The  channel  number  of  the 
maximum  count  in  the  ref- 
erence peak. 

Channels  from  the  front 
(high-energy)  edge  of  the 
reference  peak  to  the 
"third-height"  position 
on  the  front  edge. 

Channels  from  the  back 
(low-energy)  edge  of  the 
reference  peak  to  the  third- 
height  position. 

Channels  from  STOIPR  to 
the  third-height  position. 


For  MCA  operation,  the  "third-height"  position  is  de- 
fined as  STOIPR  and  hence  the  program  sets  DELTA3=0.0. 


Standard  Spectrum  Data  Set 
D1,  (Y(I) , 1*1,  10) 
FORMAT  (F6.0/10(F6.0,  2x)) 


D1 , D 
Y(I) 


Every  tenth  channel  of  the 
spectrum. 

The  number  of  counts  in  the 
Ith  channel  of  the  spectrum. 


etc 
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Peak  Position  Cards  (Use  with  NPKRD  > +1.) 

PEKPOS(I) , IFIXPK(I) 

FORMAT  (F9.1+,  II) 

PEKPOS(I)  * The  estimated  third-height 

location  of  the  Ith  peak  to 
be  considered  by- the  program, 
specified  in  channel  number. 

IFIXPK(I)  = 0 or  blank  The  program  will  adjust 

PEKPOS(I)  for  the  best  fit 
to  the  data.  (Standard 
option. ) 

> 1 PEKPOS(I)  will  be  held  fixed 

by  the  program. 

The  cards  may  be  in  any  order.  The  program  will  rear- 
range them  into  ascending  order  of  PEKPOS(I). 


Peak  Position  Termination  Card 

PEKPOS(I)  < -1.0  This  will  terminate  peak  locations 
and  the  program  will  process  a new  set  of  data, 
beginning  with  Card  1 . 
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APPENDIX  C 


RECOMMENDATIONS  FOR  INPUT 

The  following  suggestions  are  intended  to  aid  the  user  in 
setting  up  the  input  card  decks  to  ALPHAFIT. 

Backgrounds 

The  automatic  background  option  is  sometimes  capable  of 
providing  a reasonable  estimate  of  the  background,  but  it  can 
also  produce  very  unrealistic  results.  The  use  of  the  automatic 
option  is  not  recommended  if  high-quality  results  are  desired. 
Background  option  I (NBACK  > 0)  gives  the  user  the  maximum  con- 
trol over  the  background  curve.  If  enough  points  are  specified, 
the  background  can  be  tailored  to  any  specification. 

With  background  option  II  (NBACK  = -1 ) constant  backgrounds 
can  be  entered  into  the  program.  Option  I can  also  be  used  with 
NBACK  = +1  if  DB(l)  specifies  any  abscissa  in  the  middle  of  the 
spectrum  and  BACK(l)  specifies  the  background  value. 

Peak  Selection 

The  option  that  requests  the  program  to  hold  positions  of 
peaks  fixed  while  obtaining  their  areas  is  not  generally  recom- 
mended. It  is  useful  if  the  program  has  troubles  in  separating 
closely-spaced  peaks  whose  positions  are  known  in  advance.  How- 
ever, the  positions  of  peaks  are  almost  never  known  a priori. 
Even  if  their  energy  values  are  precisely  known,  experimental 
uncertainties  and  insufficient  yields  may  make  their  positions 
in  a spectrum  poorly  known. 


1 

Printout  Options 

The  input  NSWW1  = 1 produces  messages  from  the  variable- 
metric-minimization  package.  This  is  useful  for  error  tracing 
or  debugging.  Otherwise,  set  NSSW1  = 0. 

The  input  NSSW2  = 1 produces  a printout  of  the  channel 
numbers,  raw  data,  and  individual  peaks  fit  to  the  raw  data. 

It  is  sometimes  useful  if  the  user  wishes  to  make  minor  correc- 
tions to  the  final  results  from  ALPHAFIT  by  manually  summing 
the  counts  in  a peak. 


APPENDIX  D 


INTERPRETING  THE  OUTPUT 


Special  Messages 

Most  of  the  printed  output  is  self-explanatory.  In 
addition  to  the  error  messages  indicating  improper  data  (listed 
in  the  descriptions  of  the  utility  subroutines),  a few  other 
messages  can  be  produced  under  special  conditions. 

nn  Peaks  Deleted.  Some  of  the  peak  positions  that  were 
read  in  are  not  in  the  range  of  abscissa  values. 

Epsilon  Too  Small — Cannot  Converge  For  Peaks  At  xxx,  yy . 
Twenty-five  complete  iterations  have  been  made  by  the  variable- 
metric  package  in  an  attempt  to  minimize  the  function  F for  the 
peaks  located  at  the  positions  listed.  Convergence  was  not 
achieved  and  the  results  of  the  last  iteration  are  returned  to 
the  calling  programs.  Normally,  these  are  sufficiently  accurate. 
However,  this  message  is  sometimes  produced  if  the  spectrum  is 
very  complex  and  either  too  many  or  too  few  peaks  are  specified. 
One  should  revise  the  estimates  of  the  peak  positions  and  rerun 
the  program. 

Results  of  the  Peak-Fitting  Analysis 

The  final  printed  summary  of  results  should  be  examined 
for  signs  that  the  calculated  fits  to  the  data  may  not  be  re- 
liable. If  necessary,  the  input  deck  can  be  revised  and  another 
pass  made. 

Errors.  One  can  roughly  estimate  the  errors  expected  for 
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the  final  positions  and  areas  of  the  peaks  and  compare  these 
with  the  calculated  errors.  Large  discrepancies  require  atten- 
tion. The  expected  error  in  the  area  (NUMB  . COUNTS)  A is  yJ-A- 
The  expected  error  in  the  peak  position  is  roughly  FWHM/(2.5V~A) » 
where  FWHM  is  the  full-width  at  half  maximum  of  the  peak. 

Normally,  the  positions  and  areas  of  the  peaks  will  be 
determined  much  more  precisely  than  the  listed  errors  would 
indicate.  A large  discrepancy  between  estimated  and  listed 
error  is  thus  only  a sign  that  the  program  experienced  diffi- 
culties in  the  analysis,  not  that  the  answers  are  necessarily 
unreliable. 

Negative  Peaks . The  program  sometimes  produces  peaks  with 
negative  areas.  This  is  usually  a sign  that  ALPHAFIT  has  been 
asked  to  place  too  many  peaks  in  a small  region.  Some  should 
be  eliminated  and  the  data  should  be  reanalyzed.  If  the  user 
intends  to  have  all  the  peaks  he  requested,  a revised  set  of 
estimated  peak  positions  can  sometimes  eliminate  negative  peaks. 

Peaks  with  Zero  Areas  and  Errors . The  printout  may  show 
that  some  peaks  have  areas  and  errors  which  are  exactly  zero. 

This  results  from  a spurious  condition.  During  the  variable- 
metric  analysis,  these  peaks  have  been  moved  outside  the  range 
of  the  spectrum  in  the  work  area  (not  outside  the  total  spec- 
trum) and  therefore  have  zero  amplitudes.  The  errors  are  arbi- 
trarily set  to  zero.  The  cause  is  usually  the  same  as  for  neg- 
ative peaks  and  should  be  treated  accordingly. 
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Program  ALPHAFIT  Listing 
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">0  193  T=l,  JV\'< 

MJ  = MJ  4-  1 
0 = C • •? 

no  1'*]  I3  = l,Jr> 

143  n ” 0 + vimit1) 

9HVCIJI  = r,/=’L'*aT  ( J°) 

OAO<(NJ>  = . 

'■M-'K  - j-t)  ♦ M .) 

rr  (3Up?(M  J)  JoC^SC, 152 

1*5  ] M J = N J - 1 
<50  ’0  i73 

1*3  ■>"C4J)  = 3107  C-’a+il  ♦ QF~Ef.*SHIFT 
OOfMJ)  = HE(MJ) 

...TEST  r*5 3 LTM7T  OF  * C POINTS  AMP  TAKE  APOPODPIATF  ACTION... 
I'O  Te  (NJ-tl)  i5:,r,r  ,4C? 

1 7 *5  TC  C j,,-r.MO)  -T  -0  UC : 

T J = N° 
on  '0  93 

18?  ^ 5 = H ♦ j° 

20  3 r r,‘i”T'|!jr 

43  3 Tf  (mj-3)  4i5t42  3,**«C 
413  - = vj  - 1 
P-0  TO  4?C 
4?  3 •*  = 4 
4-*3  N = NJ 
533  M«  = *1  ♦ 1 

00  =13  1=1, M 
513  «(’)  = l.C 

PAI.L  LS0ODL  (*•«,  M n 
on  74*  1=1, MCMINN 
010*7(7  ) = *(1) 
on  r?c  T5si,'i 

523  0A''K?{D  = 7A0K'»rr)  + y CP+l)  *0IST  (I)  *•  IP 
TV PAC<?(  T>  ) “73  , ‘.'•+C  , 54  3 
C3G  M'^C)  = q.C 
4 «?r.  o no-riNME 

on  *0  8 33 


. . .IN"rPor>HT"  THr  n Af  K GD CUND  TM*>P!I5H  THE  “PINTS  orA0  IN..J 
?3i  j"4  = : 

J C ^ = G 


-r " J',!  = i 

1-fM  T TT  r(T  ( * ) 

( 1)  =-);\r^(  1) 

717  TCf-'T  (.J  J> -ni^”  ( irunsf!)  ) ' 1 f , 7? 0 , 75 0 

7ir  jr**i 

= V”~  ” ( 'TWA-  N ) 

(MJ4-.J71  + 1 ) ="*3r  ' ( » ■ J > 

7*0  •'O  7=1, MJ 

nrfT4.j-  1)  = m ( T ) 

7~z  sh-'i  ) = nao<  (I ) 

M='  l + J9l<-  JO? 

*r  (!KLT.7)  "O  r9  n 0 C 
J = ** 

no  ’37  1=  1 
TTV  = 313  7(7) 

773  r- niy-i- ( j*i ) ) rr>. 3f?7=, 17c 

7'n  J = ) 4-  1 

TP  (J.0E.‘J)  J=  ) - 1 

78  3 n*--MS  (0T(  J-H  -*>E  (J)  ) - r-E(  J-l) -9E ( J*i ) ) 

O*  »?  = <?!(  J)-nr(  J-l  ) ) * (nr  ( j)-C£(  J + l>> 

(n't  ( j*i ) -''T  f J-l) ) * (rr  (J*i)  -n£  ;j) ) 

3 A''*  2 (I  )=  (97<  -''r  ( J ) ) * ( 9 TV -HE ( J*i) ) /7'-1i*RH02  ( J-l) 

1 * (3Ty-r'r(J-i)  ) *•  (riy-"T  (J  + l)  )/PE‘<'!*pH02(J) 

2 «■  (3TX-?r(  J-l))  - (niX-r'fJ)  ) /9EH3*PHC2  (J  + l) 

Tc‘r»an<?C)  ._T.7.:)  *flcyp(i)=o.: 

707  * rv>7  TVJiJcr 

V'7—  £(5,^10  r'."(T),s*J02(I),I  = l,K') 
ail  e,0’MaT(2o4f.»>‘3tri7«s  |J?rr  sv  p AC'<G90UNn/ ( ?P1 2 . 3)  ) 

GO  TQ  177  0 

. . .Ca.T'JLftT-  'HP  oAO<r-  = OU‘lO  INTERVALS... 

93  7 T J = ,? 

■'*V1>  = H^Ml) 

TP  (MJ.EO.i)  90(7)  = 0197 (NCMANN)  ♦ GP7EN 
on  i2  3 1=  1 , *I?M  A ‘-I'! 

IP  (OISTd)  .IMiCJ))  TJ  = IJ  + 1 
T c (TJ.co.Mj)  m(\'J  + l)  = njGT  (NCHA'CJ)  ♦ DCT£N 
923  n/O'dC)  = 9Ar:<(;J-i) 

W°*T£((i,9?0 

Q-*3  c,r»  «S.T(1?H0  V^'A^CE  ,1'XilH  OAOKGPOUHn ) 

90  ~i*G  1=1,  MJ 

9TV  = 3Q(T*i)  - oct  t,i 

9^3  ’o—pfo,?^;)  (r»,(-),rix,^Ar<(i)) 

9K  3 c'r-'ftT(r),7,i4.,c;  .3,1-Xr6.2) 

10  3 3 

c,p 

9tn*0UTTNE  Lr'>°0,_  ( “cm , m<  ur) 

Lr>  '7  ?0J  a ” «0t.v‘>rVTM  rI~ 

<‘0110*1  /*32/  y(32).  Y(y2),  V(93),  3(10,  A(1?,1C),  WIZR  ( *C  ) , 
1 *V|  -fviy  (U  >Q) 

*.'s*'7m 
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*1- '*1 
'•’r  ‘ + M4-M 
'1  ■»«  r M ?-  ^ 
tM?1  4-M 

c-n^ATTT)  AMO  ’MVEPSICN  r r SYSTEM  Cc  MORTAL  EQUATIONS 

no  13  2 K?  = Mi , ‘"'t 
Ynnw~»f  <2)  =0.  c 
122  roirx'iu" 

nn  ?oj  <i  = i , m 

▼ ( K'l  ) 

no  200  <?  = *! t M7i 

xpt  -frof  <7 ) =t-  *•*  + Y.-*nwr J ( *?) 

>ET1=  < (<!)*’ TPM 
203  fionwE 

no  ’03  1=1, M 
nn  ’-jo  j=i,m 
*<2=I  + J*l1-1 
A (T,  J)  = y^ohe?  r-'O) 

300  cn-riN'Jt 

on  .CO  Ysi,M 

-roMrvjfi/)  »Y«  ) 

no  fjO  j 

XCO-.trp(  X2)  =-“->'1*vop'.^rc  ( x = ) 

Tr-,M  = x(<)  * 

too  •‘'n‘’Yi\'utr 

=y?o.)?°(<0) 

5 3 0 OfVIT  I'4'JE 

CALL  1ATTNV  (A  ,v*»",,lf  nETEr  '1  *13 ) 

7oo  rn,rr,i"E 

rMr 

mr'm.jTiMT 

. ..OOWrPT  n/*rA  "O  cuANv£L  HUNGERS  AMP  PREPARES  IT  FOP  THE 
•np:AHL--M  = -‘3rr;  =Ar<5Gr... 

pcf\,_  j_rr 

CC‘’*0‘l  /AGO/  <00,20  ,PcR<?0)  ,GH<20>  , <'<?0>  ,Xf»(20)  ,0®(20>  ,’(20)  , 

1 n5(23)f  co,r^tri.tsL,rp,rsc,T0,TZZ,O,  Ar,CS'',FQ,r,TP,cn,GTT,GSn, 

2 '*TLTA  , 2 E , L'  , «**,TT,L,lHCLr<C«>  •NS'^l.HSCW? 

CTMON/a:  l/*»(2i*  ",  20  ,SF  (?![?>  ,I7ST=T  (20)  , IZSTP3 <?*> , 3 < 21 0 2)  , 
l v»?i:c)  , = (2i:c) ,"(?:, 2:)  ♦ a(?d,  G<2r>,  ^t?o)f  nrLt 

? f,  f.ot«*-,  lpc(22),  lp<22a,  :pc,  :°p,  :3cf,  tic,  itc,  iof<? 

ro-'^CN  /a:c/  ctv‘A(2C),  £“:,c'7(.?j),  n"r>Tn(ri',  rTOP“-((23), 

1 V1NTI?(20>,  :-n-Y(?C),  If'TTAL(2C),  I FI  N A'_  < 2 0 > , STO?V'l  ( 210M  , 

? **A»IY{ -»i  3c>  , 1’TMEG,  ‘'NY 

rn..«QN  /AO7/  'v'  <‘'v  <-b7>  . !'•*  , \C  A'.C  , IGS,  LO 

.C'v«'*3*i  T'ITprr  (20 , 1'  ) ,ETP’’E;f  (21r  })  , n EL  T A1 , CEL  TA  2,  nEL  T A?  , 

1 'TPIOP,  prrMiY,  TO!,  I f 2 , in;,  MO,  •»',  fP.EP 
ro^HQN  /DATA/  0IrT  (213L  ) , COUIPS  < 21 J C > , ^UNIO,,  7Q0AY,  F3ET,  OFTEN. 
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1 ~-srn(-*)  , uftyoT-r  ? »mxd2 

*3.1 

no  = rjisr  (i) 

...CO MVr?*  'LL  'JUTTS  TP  CHANNEL  NMiniR... 

00  <js:  T=i,M* 

'5c"'OT0(T)  = (r'cr,'.T'i(7)  - nn)  /0  CTEN  4-  1.0  !. 

left:’  Tr’mc)  = i.' 


nc-'^TO  ( *)" ) - 0 * 1 0 . C 

">-L7  » *1  »IT  C'  rL " ; 1 ♦ 0 n_ 7 * ? C . 1 ) 

to->=  C 2 . OOFL  T \ ?) 

TO--  = { (2.34-OrL'"^)  -FLOAT  < I CR ) ) * 1 j • G 
O013c>lT  = l,,<'vn? 
ioe.i  o(T>=?.3 

on  1*1. tot 

10P."  o r VI)  = STORC-C'(T) 

TOO  = I-»r  + 1 

Tl'  = l 

iiT  = i* 

TII  = (M'lNTER  no  PEAK*  AtRZAOY  ANALY^EO)  + 1 
2313  TO-r  S III  - 1 

. ..OT'/TOO  STOMO  pc  ?rf<s  INTO  NON-OVERLAPPING  SUOGROLP  S . . . 
MS=V 
Lr=l 

TcfIi:-U’>  ?Ki,2  - ?1, 70  CG 

T10GS  = (rT»ST  ~HANNrt  NUMBER  OF  CZ?^T  PEA  < IN  SUBGROUP)  - 1 
2GC 1 TJ-SS  = orps-n (-T- ) - PTL7A?  + P.l 

7 c(  7 1C^S.  GT  . 1 ) UTSS  = 1 1 CCS  - 1 
no  >oin  jj*i,2  3 
LT=L7  +i 

TFforpp-o'  (rlZ*  IJ)-PEP.S~Q  C II I +J  J-l)  - DEL  7)  2Cir , ?0 IF , 2020 

2010  onM- i>|i|o 

2023  Y<Y  = o”RST0(T:T) 

T7C  = TOTAL  U"v"3Er  OF  CHANNELS  IN  T«E  SUBGROUP 


TT“‘  s PERSTOf LT+I0FP)  - vyx  + DEL7A1  ♦ nrLTA2  + 3.1 

...FE^JR  7 N T T t a L N'TFlv,  0A7A,  WEIGHTS,  AMO  PARAPET rRS .. . 
no  ?0?V  Tri,LT 
■»0  ??2‘  J=1,LT 
2C23  MPr:,j)Sg.c 
OH_~A  = l .3 
no  *»C?C  1 = 1,  L"p 
r“  = in pr  ♦ i 

PF-MI)  = ‘,r'3>  TO  ( IM  J - XYV  + OELTA?  + 2. 


«HCT, I) =ErCRHH  ( I J) **2 
7 c ( T HQL2 ( 7M)  ) 2:2?,2u2E,2C2? 
2C  ?■"  r C-.T 

on  "0  2:^7 

2021  RFltA  = n"L7AfcNM(Tfi) 

ao’i  oo-riMM? 

°n  ■*<} it  i = tic,.ttp 

JJ  = I ♦ li^pc  •-  4 

x r i =pa  ( jj> 


SSS55S1- 
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4 


207* 

- et_  n.\  T ( T 7 - -L  T ) * : . f ~ 

"ALL  DA  i/i  vi 

TT*  'ESs.TTI^rs  ♦! 

= i|iMr*  ir  Bfjyr  j?1  S'nGftOPP 
TMITAL  = C'T\'’TI,JG  C.J£’NEL  NvJHPEP  TN  GUt>GROUr 
I F T*l  AJ,  = CT'-'\L  NIJHOSR  FCR  "JPG^nijo 

T>n~x  (TTTMerj  =•_" 

TN'“?.L(  I7I,j-">  - Tl?c' 

Te,T\’AL  ( TTTMrr ) = TlG"r  ♦ IT"  - 1 
Vs  = IT z - l: 

...TALC.  ~:'~2  pG  TM  PCSIT:CmS  ArjO  COM^Tf  MO*  HA!  I7FD  I M "?  I V • PEAKS 
"0  20V2  1-1*1.“ 
tm  = T?rr  + T 

r^GTO  ( J14)  = (v;  ( t)  ♦xxy-OELT  A2-2.  C)  *QrTEN  ♦ 00 

CT«->t)w»(  T -J)  -m  ( r,  y ) 

(IH>  s?r;»T  ( - -> c ( (2  .OSTORHH  <IH)  *f^eS7)  /Q?)  )*OFTEM 
Tr(*,(T))  ?li",  711-T  f 2115 
211  * -517(!M)  = :.c 
GO  *0  2120 

21  lr*  (IH)  = ^fl=‘.(T)  -/A  (II 

212:  GC  •’•Ji*:  jsi.-avs? 

?0tO  ■’(1,1'*’)  = 7 ( J , I 4 ) *•  £ ( I ) 

■;  *>  2 rryriNOE 

T r (IT. OS. 2")  ••»2T“r  < e , 1 C C i ) (PEPSTO(T) ,I=1,LT) 
ttt  •-  j xi  ♦ LT 
on  TO  2:3  0 
300  3 D~'RM 


10G1  rO->  J£T(  Oimj  r n G I L O N TOO  HI  ALL  - CANNOT  CONVERGE  FOR  DEA<S  AT  ,5F10 
1.2'(15c1j.2») 

HP39'ITINE  nsnr-> 

OO'HON  /A  3 C / H»(2>:  , 20)  ,PFP(2t.)  ,GR(20  ,F(20'  ,XP(2G)  ,0=(23)  ,T(23)  , 

1 GP  ( 20  > , - , G"  , EL  . SL  , F = , r,«-e  , T 1 , 7 ,0  , A , OSS  , F 2 , GTP  ,F« , GTT  , 959  , 

2 L TA  ,F,'V4?,IT,  L,  THOLE  (22)  , NGSW1 * f'SSW2 

jn  = N - 1 

n0  100  1=1, JO 
K = I ♦ 1 
10  4 30  J*K,M 

tc  (pep (i > -osn ( j ^ ) ic:,i'c,s: 

SO  FAV:  = °ER(I) 

rc-d)  = pcs(j) 

ocr^(J)  = SAVE 
GA = MH(T,f) 

HH(T,I)  = J»(J,J) 

MH(J,J)  = 

T«*AV  = iHOLOd) 

TM->i_n(i)  = iholkj) 
iwfLO(J)  = TSAV 
133  r n*** i *;m E 


»r“  !RM 
“NT 


■ _e<r  QUA^^,, 
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‘?l|V>'Vr  jvjr  c-»,(LL'rT  ,r.c  f rr  f err, Ml ) 

, . . *"T^  J3  *-C*  i.IV-fiP  /fit  AT?(/1S  AN0  C ?,  TA  T NS  PELA’IVE  »*»K  HEIGHTS. 
a*'L  LpP 


nT'-NSITI  rT?  C>  3 > , Op  (2  • > 

~C  • «ON/ATl/^f  ?i:  ?,  ?:•>  rcir?)  ,IZST^T  (20  ,*?STPP  (23)  ,?(  2i  3?)  , 

1 *(210 3)  ,<-.,<?m)  , m?o,  r; < ?c > • ”?.<2-j>  » rs.» 

2 n-P,  F,  pnre,T.  L3!T(?2),  LP(22),  IP*?,  IP?,  In?E,  lie,  ITC,  IOr? 

/Tin/  ot*t (?i: ’ ) , roller? (213 o , ^>ur!in,  today,  =-»eo,  octen, 

1 ,l'3rT?N»  MOH.VJM,  II,  IJ,  IS°SCT,  IPMT,  lE'-'R,  MAXPEF,  ucyP2 

L**L17* 

TFf^l.-n.l)  pn-S*  = l.CEiO 
fcC06  0 0 11  L=1 , l.T 

L °(L ) = 3E'(L) 

L5r(L)  = (P-f  (L>  -'LO^ILMD  > ) *13.  i 
91  CALL  SHT--(L) 

2Ce»Tp(Ml-3)2?2,21‘>,,2? 

O'JMMY  0A Lu  tO  c,.IN  to  pO?H  THE  A??AY  F 
210  9 ALL  r,,N(  L”) 
no  v T=IiC, 

220  W(-)  = l.:/F-!5T(A9F(p(I))+i;.) 

22?  nO  »i:  >1=1, LT 

MM  z M + I O F - 

"9  210  L =1,H 

ll  - l io~r 
* csl  ) = j . : 

no  203.  Islin,  1*0 

TF  ( (I.LT.t-'-TvT  fi_L)  ) .0?.  n.5v,I7CTOr(LL>  )1  00  TO  233 

T*  ( (t.lt.t7ot->t  ( jh)  ) ,oo,  (i.gt.:zsto“(ym>  n go  to  2 j 0 

TL  = I - I79T?T(lL)  +1 
™ = I -I2FOT  ( 1*0  ♦ l 

CCSL)  = C(M,L)  ♦ H(I)*W(I)*Z(IL»LL)47(IH,mD 
20  3 POHTINME 
21C  0 (L,M>  = 0 CSL) 

00  ?01  H=1,LT 
MM  — ♦ TC-e 

a < »>  = a.? 

00  231  1=110, t-o 

Ip  ( (I.  LT.I^S^nr  (m.m  ) ,0R.  (I.r,T.I7ST0r(MM)  ))  RO  TO  201 

TM  = T - T7G“^~  «•  1 

ACM  = ?.(M>  * U I)*H(l!»y(I)‘Z(IM,l^) 

2C1  OPMtTnUE 


CALL  M A *T  v|  V ( 0 , L T , A ,1  ,OFTrp^t?]) 

r ALL  nF?TVD(L*toF3) 

no  160  0 L * 1 . L T 
190  3 OM)  = 0 ( L ) 

CALL  =XCJ('_T) 

’?mi-t)  2r!»iUii,?oe 

ltll  r,T  = -/  CTTC-I10-LT) 
no  1306  L = 1 , L T 
»-•»  = ”7*3  (L  , L) 

1C 0 6 p^n(L)  = SOP*  C\9G(pE?> ) 


SSSSSS^ 


I 


▼ r ( r ? , LT  . rr,"CT)  '■'’ETT  = 

3 C-i13 1 ij 

T ,p 

gm^QijT  T»;r  nrrr  ,r'(l.(_TT>  rrP)  i 

• • A LC  t)L  1 Te  = ->^rI7iTIv/E5  OF  FUNCTION  F(FU‘l)  WIT^  CHIM^F 
VI  f'C'5I7ln*' 

->'*L  LP' 

nT  i-  *!«;:  tj  B-’  (•>  j) 

• r r*^0*l/«3  l/’(  V ? ’’j  71  > , Fe  r?  - C2>  ,I7S‘r2T  (EO)  ,T7STP»  (20  , ?.(  2102)  , 

i y c 21c  o ,wr»i  n> . * (?i:  j>  ,c(23 ,20  , m?  t>,  ^ ( 2 c ) • e2R(70),  nrL, 

*»  F,  F’*:'**,  I -F(<2),  IP(72),  I'7  ** , IPR,  IPr",  112,  ITC,  I2Fr 

BT'-JON  /TIT!/  T[~-~  (21  w-.  ) , F OI'NTF  ( 7 10 !.  ) , P'^'in,  TCO'’^,  roE'),  3CTFM, 

1 ,'-»FTE‘J,  \'C-U'<'I,  :i,  IJ,  ISPEOT,  ICHJ,  IE^p,  Mf  VPEF^ayP?  •. 

L“-LLTT 
nn  -331=1,17 
u.  : L * TOrr 
no  P'. 

gt,  <*«;(-)  = 7(1, L'O 
T«f'7i  = its-b-(LL) 

T V«/2  = I7«;to=(lL) 

LO‘,=  L°(L) 

LO**-  = l°F(L) 

Y=p-b<L)OpL 
L p (21 ) = V 

L °r  ( 2 1)  = (Y-FLT  5T ( L° ( 2 1 ) ) ) *13 . 3 
v=->  = ?(L)-1EL 
L c ( 22 ) = V 

1.  P"  ( 22)  = (V-C1.0.1T  U.P(22>  ) >*la  .3 
L °(L) =L  °( 21) 

L0'<L)  = L3r(?i  ) 

^ ILL  SMTeT(L) 

COM.  r! ) N ( L T ) 

Fr=  = 

l°a)=LD(22) 

L Pr  ( L ) =LDC'  ( 72  ) 

C li  t.  Fl,?rT  (L) 

* ILL  piJN(LT) 

r.  (._>  :|ro.r,  /(?.]*npL) 
n n '» r,  j =i  f m * v o 7 
= t^(T) 

I’O’gTCLL)  =17  171. 

I7~.TC°('-_)  = TF1V2 

L°a)=LC)P 

Lor(L)*LP»E 

np 

P«r^OUT?Mr  euTr-(L">) 

. ,,"0'.'T 70'. 7 :V‘rE*,PCLiTICt  PF  ”c,  PEAK  POd  SHIFTING  Pr>1CTI0NS 

nr  CHANNEL? 
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L°  c 

~ 0 ION  / \Z  i/-'  ( '‘IT,  '»•;  ) , c ~ < Pl„2)  , I7ST-?”  ( 2 , ) ,*Z0T'P  (C)  , P ( 21 . 

i v(  zv:  :> , w<  ?i : I ) , r(  ?i : :>  ,c  <?c  » 2: ) , m?:),  0(2  :>,  r??.(2;), 
? *>-o,  =•,  FT"**',  |_  “ r ( 2 ? ) , L P ( ? 2 ) » IV,  T^,  I3'517,  11:,  I’( 

como*'  nw  "r cv..-:)  ♦ covers ( ? 1 ;» r ) , pyn,  cvf 

i H-Hi  |>|,  II,  IJ,  IS-ECT,  Ir?iT,  Troo.MftxprP^SVPS 


L=L~> 

LL  = L ♦ I0r- 
no  13  2 JsltM/ix3? 
r.  I -L  3 (L  ) 


n r>~ i_oc  ( (_) 

IJ='1*0?/13  .3 

n q=  T 35 


0 O-T  o.'jr 


1C  3 "*  ( J , L L ) = ?.C 

01=  -»c  +10/ 1C  . C 

( 3l-'J)  ?0C,1J*° 

0*P  VriNP  P'-<Ir“L  L L£S0  ThAH  3 
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